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Abstract 

We introduce a simple cellular-network model to explain the capacity of the plants 
as memory devices. Following earlier observations (Bose [1] and others), we regard 
the plant as a network in which each of the elements (plant cells) are connected via 
negative (inhibitory) interactions. To investigate the performance of the network, we 
construct a model following that of Hopfield, whose energy function possesses both 
Hebbian spin glass and anti-ferromagnetic terms. With the assistance of the replica 
method, we find that the memory state of the network decreases enormously due 
to the effect of the anti-ferromagnetic order induced by the inhibitory connections. 
We conclude that the ability of the plant as a memory device is rather weak. 
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1 Introduction 



Since the pioneering work by an Indian scientist J. C. Bose [1], plants have 
been regarded as networks which are capable of intelligent responses to en- 
vironmental stimuli. For example, the dodder coil, which is a plastic plant, 
explores new host tree within hours after their initial touch contact [2]. This 
sort of behavior might be regarded as plant intelligence. If that is the case, does 
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the plants compute, learn or memorize various spacial and temporal patterns 
in different environments as computer or our brain does ? 

Recently, Peak et. al. [4] pointed out that the plants may regulate their uptake 
and loss of gases by distributed computation. As well known, the ability of 
neural networks, which is a mathematical model of brain, is also based on 
parallel and distributed computation. Therefore, similarities between neural 
network model of brains and the plant network should be discussed. Although 
the behavior of the dodder coil we mentioned above is due to emergence of 
the intelligence as a macroscopic function, it is important for us to investigate 
its microscopic reason. 

Almost eighty years ago, Bose [1] detected electrical signaling between plants 
cells. Since his experiments, many examples of cross-talk, namely, the bio- 
chemical signaling pathways in plants have been found. Especially, a Boolean 
representation of the networks of signaling pathways is possible in terms of 
logical gates like AND, OR and XOR etc. These Boolean descriptions make it 
possible to draw analogies between plant networks and neural network models. 

Recently, Bruggemann et. al. [5] found that the plant vacuolar membrance 
current- volt age characteristics is quite nonlinear and almost equivalent to that 
of a Zenner diode. Inspired by such observations, Chakrabarti and Dutta [6] 
identified the cell membranes as the two state neurons of the plants and utilized 
such a threshold behavior of the plant cell membranes to develop or model 
gates for performing simple logical operations. They found that the plant 
network connections are all positive {excitatory) or all negative (inhibitory), 
compared to the randomly positive-negative distributed synaptic connections 
in real brains. As a result, the plant network does not involve any frustrations 
in their computational capabilities. Hence, although the logical gates for com- 
putations could be achieved by such networks of plant cells, for intelligence of 
such networks they must must also possess some memory capacity (in order 
to caompare and optimise). 

With this fact in mind, in this paper, we investigate the equilibrium prop- 
erties of the Hopfield model in which both ferromagnetic retrieval and anti- 
ferromagnetic terms co-exsist. The strength of the anti-ferromagnetic order is 
controlled by a single parameter A. Within the replica symmetric calculation, 
we obtain phase diagrams of the system. The A-dependence of the optimal 
loading rate ac at zero temperature (T = 0) is discussed. 

This paper is organized as follows. In the next section, we introduce several 

experiments and observations about the current-voltage characteristics of the 
plant cell membrance. In section 3, we model the plant with such properties by 
using a Hopfield-like network model in which both ferro-retrieval and anti-ferro 
terms exist. In section 4, we analyze the model with assistance with replica 
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method. In this section, we investigate to what extent the ferro-magnctic re- 
trieval order remain against the anti-ferro disturbance. We also investigate the 
result of the ferromagnetic disturbance. In the final section is summary. 



2 The I-V characteristics of cell membrances 

In this section, we briefly mention several results concerning properties of 
the plant units, namely, current (I) - voltage (V) characteristics of their cell 
membrane. In Fig. 1, we show the typical non- liner I-V characteristics of cell 
membranes for the logical gates. 



Fig. 1. (a) The non- linear I-V characteristics of cell membranes (of. [5]). (b) Its 
Zener diode like representation with threshold voltages denoted by vt- 

From this figure, we find that the I-V characteristics is equivalent to that of 
Zenner diode. From the view point of input-output logical units like percep- 
trons for neural networks, the output of ith unit Oi is given by [6] 



where the connections Wij are all positive or all negative. In the Hopfield 
model they are ± randomly distributed weight matrix as given by the Hebb 
rule. From these experimental results and simple observations, we now have a 
natural question: could the plants act as memory devices as a real brain does ? 
Obviously, in the above definition of a single unit of the plant network, there is 
no frustration as in the animal brains. Thus, a scope of this paper is to make 
clear this problem, that is to say, to what extent, this kind of limitations 
for sign of the weight matrix influences the ability of pattern retrieving as 
associative memories. 
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For this purpose, we introduce a simplest plant intelligence model based on 
a Hopficld-Iikc model in which ferromagnetic retrieval and anti-ferromagnetic 
ordered phases co-exists. In the next section, we explain its details. 



3 The plant intelligence model 



We start from the following Hamiltonian : 




(2) 



'Haf = -TT^SiSj, Hfr = -jr^^^i^jSiSj 



(3) 




where = (^f , ■ ■ ■ , is fith embedded pattern and S = {Si, ■ ■ ■ , Sjy) means 
neuronal states. A single parameter A determines the strength of the anti- 
ferromagnetic order, that is to say, in the limit of A — > oo, the system is 
completely determined by Haf- On the other hand, in the limit of A — > 0, the 
system becomes identical to the conventional Hopfield model. The purpose of 
this paper is to investigate the A-dependence of the system, namely, to study 
the A-dependence of the optimal loading rate etc (A) at T = by using the 
technique of statistical mechanics for spin glasses. 



4 Replica symmetric analysis 

In order to evaluate macroscopic properties of the system, we first evaluate 
the averaged free energy : 



where <S • • • » means the quenched average over the p — Na patterns. To 
carry out this average and spin trace, we use the replica method [9,10]. 

After standard algebra [9,10], we obtain the pattern-averaged replicated par- 
tition function as follows. 




< log Z >^ = -Clog tr{s}e '^^ >^= lim 



n 



(4) 
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« ».= n / , ' \ ; " (5) 



oo +ioo 

dr, 



X J dq^p J -^eM-Nf{m,q,M,r)]. (6) 



-oo 



By assuming the replica symmetric ansatz, namely, — M, rUa — m, — 
Q, fap — f-i we obtain the free energy density per replica number n as follows. 



f(vi.q.M.r) 13 2 1^^ 2 , N 

= 0^ +^3-1-? 

n III 

H-|{logll-/3(l-,)]-^-^) 

00 

- log ^ Dz log 2 cosh P{Xm + y/arz + M) (7) 

—00 

where we defined Dz = dze~^'^ ^'^ / \/2t:. We should keep in mind that physical 
meanings of m and M are magnetization of the system, overlap between the 
neuronal state S and a specific recalling pattern among aN embedded 
patterns, respectively, q means spin glass order parameters. 

In the next section, we evaluate the saddle point of this free energy density 
/ and draw phase diagrams to specify the pattern retrieval properties of the 
system. 

5 Phase diagrams 



In this section, we investigate the phase diagram of the system by solving the 
saddle point equations. 

Saddle point equations: By taking the derivatives of / with respect to M, m, r 
and g, we obtain the saddle point equations. 

00 

M 



= j Dz tanh - \)M + z^] = -m (8) 

—00 
00 

q^ j Dzts^uh^ (3[{l-\)M + z^/^\ (9) 



—00 



(10) 
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We solve the equations numerically to obtain the phase diagram. 

T = noise-less limit: We first investigate the T = limit. In this limit, 
obviously, — > 1. After some algebra, we find that the optimal loading rate 
etc is determined by the point at which the solution of the following equation 
with respect to y vanishes. 

{l-X){l-2H{y)} (11) 
Dz. In Figure 2, we plot the optimal 
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Fig. 2. The optimal loading rate Uc as a function of A. ac{\) decreases monotonically. 

loading rate etc as a function of A. From this figure, we see that the optimal 
loading capacity q;c(A) monotonically decreases to zero as ~ (1 — A)^. This 
means that the fcrro-magnetic retrieval order was destroyed by adding the 
anti-ferromagnetic term to the Hamiltonian. Thus, we conclude that if the 
weight matrix of the networks is all positive, the plant intelligence model does 
not act as a memory device. 

Spin Glass-Para phase boundary: Before we solve the saddle point equations 
for T 7^ 0, it should be important to determine the phase boundary between 
the spin glass and para-magnetic phases. The phase transition between these 
two phases is first order, by expanding the saddle point equations around 
M = g = 0, we obtain 



(12) 

(1 — X)I3 and T = we obtain 



rsG = (l-A)(l + v^) (13) 
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y{^+\-{l-X)e 



TT 



where H{x) is defined by H{x) 




By solving this equation with scaUng /5 
the phase boundary line : 
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Phase diagrams for T ^ 0: Here, we investigate the phase diagram for T ^ 
by solving the saddle point equations (8), (9) and (10) numerically. From 




Fig. 3. The phase diagram of the system. Ferro-retrieval phase shrinks to zero 
as A increases to 1. The para-spinglass boundary is analytically obtained as 
Tc = (1 — A)(l + ^/a) and is independent of A. 

Fig. 4. The overlap M as a function of a for the case of A = 0.2 at temperatures 
T = 0.05,0.1 and T = 0.3. 

this figure, we find that the ferro magnetic retrieval phase shrinks to zero as 
the anti-ferro magnetic order increases, namely, A — > 1. The behavior of the 
overlap M as a function of a is shown in the right panel in Figure 3. The 
overlap M becomes zero discontinuously at a = Oc- 




Fig. 5. The phase diagram of the system for negative A = —0.2. 
Fig. 6. The optimal loading late Uc as a function of A (< 0) at T = 0. 



Negative X case: Here, we consider the case of negative A. From the Hamil- 
tonian, we find 

'H = -lT.t^^^^s.s,-^Y.s^Sr, A' = -A(>o). (i4) 

When A increases, the system changes to the pure ferromagnet. Let us think 
about the limit of A — * — oo in the saddle point equation (8). Then, the 
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term (1 — \)M appearing in the argument of tanh[/9(- • •)] becomes dominant, 
namely, (1 — A)M ^ z^/ar even if the loading rate a is large. Consequently, 
equation (8) leads to 



oo 



Dztanh[/5(1 - A)M] = tanh[/5(l - A)M] 



(15) 



— oo 



If the factor 1 — Ais large enough, the term tanh[/5(l — A)M] becomes sgn[/5M] 
and saddle point equation (8) leads to M = sgn[/3M]. Apparently, this equa- 
tion has always a positive solution even if the temperature T = (3~^ is large. 
In this sense, the factor (1 — A) has a meaning of temperature re-scaling. It is 
also possible for us to understand this result from the different point of view. 
In the saddle point equation : 



the second term appearing in the argument of tanh, z^Jar means cross-talk 
noise from the other patterns (/j, = 2, - ■ ■ ,p) and obeys Gaussian distri- 



inant, the system cannot retrieve the embedded pattern Usually, r in the 
second term grows rapidly as T increases. And obviously, if a increases, the 
noise term Zy/ar also increases. As a result, the signal part (1 — A)M becomes 
relatively small and the system moves from the retrieval phase to the spin 
glass phase. However, if A is negative large, the signal part is dominant and 
the noise part becomes vanishingly small. This is an intuitive reason why the 
optimal loading rate increases for negative A. In Figure 6, we plot the optimal 
loading late function of A (< 0) at T = 0. As we mentioned already, 

the optimal loading rate ccc monotonically increases as A goes to — oo. 



6 Summciry 

With the identification [6] of plant cell membranes as neurons, and demon- 
strating [6] the possibilty of logical gate operations in these cellular networks, 
the discussion on intelligence of a plant could not be complete without the 
demonstration of memory in such, either excitorily or inhibitorily connected, 
cellular networks. In this paper, we introduce a simple model based on a 
Hopfiled-like network for explaining the memory capacity of the plant cell 
network. Following the experimental observations [5], we construct a Hopfield 
model in which ferromagnetic retrieval and anti-ferromagnetic terms co-exist. 



oo 




(16) 



— oo 
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The strength of disturbance of pattern retrieval by the anti-ferro order is con- 
trolled by a single parameter: one can easily see that the internal fields Haf and 
hpR contributed by Ti-AF and TIfr respectively are of order A and \/pN /N = 
^/a in (3). We find that the anti- ferromagnetic order prevents the system from 
recalling a pattern. However, even when A is greater than ^/a, the network 
still possess considerable memory capacity (see Fig. 2). This result means that 
the ability of the plant as a memory device is rather weak, if we set all weight 
connections to positive values. Our analysis here has been for fully connected 
networks. For real plants, the cell membranes should be located on a finite 
dimensional lattice [11] or on a scale-free network [12]. Investigations for these 
situations will be made in future. 
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